-IDENTIFICATION DIVISION
-PROGRAM-ID ChargeInEMfield
-VERSION 20240116
-COMMENT A mass m with charge q>0 is located in a static electro-magnetic field [0,0,E], [0,-B,0].
-COMMENT Find the trajetory
-COMMENT Differential equation (Lorentz force): x'' = omega*z', z'' = a*omega^2 - omega*x'
-COMMENT good start settings are omega=0,92, a*omega^2=0,05
-
-ENVIRONMENT DIVISION
-ENGINE Anabrid-THAT
-TIMEBASE 1ms
-REQUIRES INTEGRATOR 4
-REQUIRES INVERTER 1
-REQUIRES COEFFICIENT 3
-
-DATA DIVISION
-OUTPUT OUTPUT.X x
-OUTPUT OUTPUT.Z z
-
-PROGRAM DIVISION
-x'' -> INTEGRATOR.1 -> -x'
--x' -> INTEGRATOR.2 -> x
--x' -> COEFFICIENT.6 -> -omega*x'
-+1 -> COEFFICIENT.7 -> a*omega^2
-a*omega^2, -omega*x' -> INTEGRATOR.4 -> -z' # input is z''
--z' -> INTEGRATOR.5 -> z
--z' -> INVERTER.1 -> z'
-z' -> COEFFICIENT.5 -> omega*z' = x''
-
-OPERATION DIVISION
-MODE REPEAT
-OP-TIME 0,21ms